Scattering theory of current-induced forces in mesoscopic systems 
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We develop a scattering theory of current-induced forces exerted by the conduction electrons of 
a general mesoscopic conductor on slow "mechanical' degrees of freedom. Our theory describes the 
current- induced forces both in and out of equilibrium in terms of the scattering matrix of the phase- 
coherent conductor. Under general nonequilibrium conditions, the resulting mechanical Langevin 
dynamics is subject to both non-conservative and velocity-dependent Lorentz-like forces, in addition 
to (possibly negative) friction. We illustrate our results with a two-mode model inspired by hydrogen 
molecules in a break junction which exhibits limit-cycle dynamics of the mechanical modes. 
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Introduction. — Current-induced forces play a central 
role in many contexts of modern condensed matter 
physics such as molecular electronics [IJ, nanoelectrome- 
chanical systems (NEMS) \2\, spintronics or back- 
action forces in quantum measurements |1HB]. In its 
simplest incarnation, the current flow through, say, a 
carbon nanotube suspended over a gate electrode gov- 
erns the electronic occupation of and thus the capaci- 
tive force exerted on the flexible nanotube. This force, 
including the dissipation of the mechanical motion re- 
sulting from retardation effects, was recently observed 
in seminal experiments [3 |S]. The nature of the inter- 
action of conduction electrons with collective degrees of 
freedom, such as vibrational modes in NEMS or the mag- 
netization in spintronics devices, depends sensitively on 
the relevant time scales. In the context of NEMS and 
molecular electronics, much work has focused on the limit 
in which the vibrational degree of freedom is fast com- 
pared to typical time scales of the conduction electrons. 
In this regime, the electron- vibron interaction causes vi- 
brational sidebands in the current-voltage characteristics 
[Ij, Franck-Condon blockade of transport |j9j, or current- 
induced heating of the vibrational mode [T]. 

Here, we consider the opposite regime in which the col- 
lective mode is slow compared to electronic time scales. 
This limit applies for instance to the flexural modes of 
suspended graphene membranes and carbon nanotubes 
[7], the current-induced electromigration of nano-scale 
islands on Ag surfaces PD], the cooling and ampliflca- 
tion of mechanical motion by the mesoscopic backaction 
of electrons [5, 151 , certain molecular switches [P, TP , spin 
transfer torques and current-induced domain wall motion 
in ferromagnets |3] , as well as NEMS near continuous me- 
chanical instabilities jl2|. In this nonequilibrium Born- 
Oppenheimer (NEBO) limit, it is typically appropriate 
to describe the collective mode cLS £L slow classical de- 
gree of freedom exhibiting Langevin dynamics. The cen- 
tral subject of this paper are the current-induced forces 
governing this Langevin process. Given the wide appli- 
cability and physical transparency of the scattering ap- 



proach to quantum transport in nanoscale systems |13j, 
it seems highly desirable to develop a scattering theory of 
current-induced forces for a general phase-coherent meso- 
scopic conductor. Such a theory is the main result of the 
present paper. 

When developing a general theory of these current- 
induced forces out of equilibrium, it is imperative to con- 
sider more than one collective coordinate [T^[T5] . Indeed, 
this introduces several qualitatively new and remarkable 
features into the Langevin equation 

or r 

^ I/' 

for the collective coordinates X^. The left-hand side 
describes the purely elastic dynamics and the right- 
hand side collects the current-induced forces, including 
the fluctuating Langevin force ^^(X.). In the presence 
of an applied bias, the current-induced average force 
^^^^(X) will generally be non-conservative. Moreover, the 
velocity-dependent force does not only contain a fric- 
tional contribution, corresponding to the symmetric part 
of the tensor ^i,^> (X) , but also a Lorentz force (or Berry- 
phase contribution [16_), arising from its antisymmetric 
part. We will see below that these features can cause 
unconventional collective dynamics and may provide the 
operating principle of a nano-scale motor. 

Model. — Our starting point is a simple yet very gen- 
eral model of a mesoscopic conductor [13j . The "quantum 
dot" is described in terms of M electronic orbitals with 
fermion operators d„ obeying the Hamiltonian Hu ~ 
'Ylinn' ^rX^D)rm'dn' ■ The conduction electrons are cou- 
pled to a collection of TV "slow" collective degrees of free- 
dom X(i) — . . . , Xtq\ which may, for instance, rep- 
resent vibrational modes or atomic locations. We take 
this coupling to be linear, hjj = /io+X]!^i ^i^-^^j charac- 
terized by the coupling matrices A^. As usual in applica- 
tions of scattering theory to mesoscopic transport, we as- 
sume that the electrons are otherwise non-interacting. In 
addition, the total Hamiltonian H = Hd+Hl+Ht+Hx 
models the leads (labeled by a and with chemical poten- 
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tials fj,a) as free Fermi systems, Hl = J^rji^v ~ A'a)c3,c^ 
with composite index 77 = (k, a) , in terms of the lead 
fermions with Fermi-Dirac distribution fai^ri)- It also 
includes hybridization between leads and quantum dot 
through Ht = J2rj,n clWjjndn+ii-C. with coupling matrix 
W, and describes the uncoupled mode dynamics through 
the Hamiltonian Hx = E,. P'i/'^M^ + C/(X). 

N on- equilibrium Born-Oppenheimer approximation. — 
Generalizing previous formulations for iV = 1 [TT], we 
start by deriving Green's function (GF) expressions for 
the current-induced forces in Eq. ([ij within the Keldysh 
formalism, based on the Heisenberg equations of motion 

dU 

M,X, + Q^=-Y1 di{K)nn'dn-- (2) 
^ n.n' 

(An alternative derivation proceeds in terms of Keldysh 
functional integrals, following e.g. Ref. |T7].) In the 
NEBO approach, we average the force operator appear- 
ing on the right-hand side over time intervals long com- 
pared to the electronic time scales but short compared 
to time scales of the collective dynamics. In addition 
to a fluctuating force discussed below, this yields a 
time-averaged force tr[jAi,C/<(i,t)]. Here, the trace acts 
in dot level space, and the lesser GF is — 
i{d\,{t')d„{t)). Note that this electronic GF should be 
evaluated for a given time dependence X(t) of the collec- 
tive coordinates. For slow collective dynamics, to linear 
order in the velocities X^, we find (see below) 

g< ~ G< + ^ {d^G<A ■ XG^ - G<A • Xa^G^^ 

+d,G^A ■ XG< - G^A • Xa,G<) (3) 

in terms of the frozen (adiabatic) GF G(e, X) obtained 
for fixed configuration X. Inserting Eq. ([3]) into the time- 
averaged force given above, we obtain GF expressions for 
the current-induced forces in Eq. ([ij . The mean force 

Fu{y^) = - j |^tr[G<(6,X)A,] (4) 

follows from the strictly adiabatic (frozen) limit ~ 
G^. The first non-adiabatic correction of Q"^ gives rise 
to the velocity-dependent term in Eq. ([T|). Decomposing 
7 = 7** -|- 7° into symmetric and antisymmetric contribu- 
tions, we obtain for the frictional contribution 

7:.'(X)=y"gtr{A,G<A,,a,G>}^ (5) 

and for the effective orbital magnetic field 

7.V(X) = -y" ^tr{A.G<A,-a,(G-^-fG^)}^. (6) 

Here, we employed G> = G< -I- G^ - G^ with G^ = 
{G^)^ and the shorthand {A^y^}^/^ = [A^y- ± A^,^)/2. 



The fluctuating force in Eq. ([IJ can be obtained 
from the (thermal and nonequilibrium) fluctuations of 
the force operator d'l{K^)nn'dn' in Eq. d2j). Ex- 

ploiting Wick's theorem, one readily flnds that the 
fluctuating force is a Gaussian white noise variable, 
{^v{t)iv'{t')) = Di,y'{'K.)5{t — t'), with symmetric vari- 
ance matrix 

D,,, (X) = y J tr { A,G< A,,G> }^ . (7) 

The correlations of S,^{t) are local in time since in the 
NEBO limit, the electronic fluctuations involve much 
shorter time scales than the collective dynamics. For 
N = M = 1, our expressions recover previous results 
[TT| . We emphasize that all current-induced forces in Eq. 
([ij are nontrivial functions of the collective coordinates 
X as the latter enter into the frozen GF G(e,X). 

We briefly sketch the derivation of Eq. ([s]). Follow- 
ing Refs. [Ill HH], we start from the Dyson equation 
g^{idf - - S^) = 1 for the retarded dot GF. The 
coupling between dot and leads is accounted for by the 
self-energy = —i Fq in terms of the hybridization 
matrices Fq, = irvW^IiaW (j/ is a lead density of states 
and He a projection operator to states in lead a). The 
Dyson equation can now be solved to linear order in X 
by passing to the Wigner representation. Relating Q"^ to 
Q^- using the Langreth rules yields Eq. ([3|. 

Scattering theory of current-induced forces. — The scat- 
tering approach to quantum transport in mesoscopic sys- 
tems, as epitomized by the Landauer-Biittiker formula 
for the conductance, is formulated in terms of the S- 
matrix describing the single-particle scattering by the 
nanostructure of electrons incident from one lead. For 
the model considered here, the frozen ^-matrix is readily 
related to the frozen retarded GF G^{e,lC) through [l3] 

^(e, X) = 1 - 2mvWG^{€, iqw^ . (8) 

We assume the wide-band limit in which v and W are 
energy independent. The average force F^{X.) in Eq. CI 
can be written in terms of the frozen S'-matrix [Eq. (IsF] 
by expressing the lesser GF G^ = G^S^G'^ in terms 
of the self-energy S<(e) = 2i^^ faXa HH]. Using the 
identity 2TTivW'^W ^ (G^)~i - [G"^)'^, we then flnd 




for the mean force expressed through the S'-matrix. Here 
the trace "Tr" acts in lead-channel space. By a similar 
set of steps, we can also express the variance [Eq. ^] of 
the fluctuating force in terms of the frozen 5-matrix, 

D,,,iX) ^Y.j (10) 
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By going to a basis in which D is diagonal and using 
IIq = n^, we find that Z? is a positive definite matrix. 

To express the velocity-dependent forces in terms of 
the scattering matrix in general nonequilibrium situa- 
tions, we need to go beyond the frozen scattering matrix 
S. The relevant concepts have been developed in the 
context of the scattering theory of quantum pumps |18| . 
For adiabatic parameter variations, the full iS-matrix of 
mesoscopic conductors can be expressed in Wigner repre- 
sentation as S{e,t) = 1 — 2iTivWQ^\e,t)W\ which nat- 
urally extends Eq. ([s]). Expanding S to linear order in 
the velocities X of the adiabatic variables, Moskalets and 
Biittiker [TH] introduced an A-matrix through S{e, t) ~ 
5(e, X(t)) + Mt)A,{e, X(i)), where 

A^(e, X) = T:vW{d,G^KG'^ - G^Kd,G^)W'^ (11) 

in terms of the frozen GF. It is now tedious but straight- 
forward to write the damping coefficients in Eq. Q as 



E 



dS"^ dS 



(12) 



Similarly, the effective magnetic field in Eq. ([6| becomes 



E 



de 



2m' 



dA, dAl 



S 



(13) 

Several comments on these results are now in order. 

(i) The current-induced mean force F^{X) in Eq. ^ 
is conservative (the gradient of a potential) if its "curl" 



dF, OF,, 



dX,, dx. 



E 



frj 9S \ /i.N 



vanishes. This is indeed the case in thermal equilibrium. 



where Eq. ( 14 1 can be turned into a trace over a commu- 



tator of finite-dimensional matrices due to the relations 
fa = /, J2a^a = 1, and unitarity S ^ 1 implying 
d{S^ S) / dXi, — 0. In general, however, the mean force 
will be non-conservative in out-of-equilibrium situations. 
Equation (|9| for F(X) provides an S'-matrix represen- 
tation for the "water wheel" [Hj or "electron wind" [T51 
force discussed recently. 

(ii) The second contribution to the damping matrix 
(12 1 is a pure nonequilibrium term. In fact, it can be 



shown to vanish in equilibrum using unitarity of as 
well as S, where the latter implies [TH' AS^ + SA^ = 
{i/2)[dxSdeS'' - d,Sd:>^S'']. The first term in Eq. (Ill 



is a close analog of the S'-matrix expression for Gilbert 
damping in ferromagnets |19| . From this perspective, our 
expressions for 7* can be viewed as mechanical analogs 
of Gilbert damping, generalized to arbitrary nonequilib- 
rium situations. While the first term has strictly pos- 
itive eigenvalues, the second term can have either sign 



and even cause negative eigenvalues of the full dissipa- 
tion matrix 7^*. (For a simple model exhibiting negative 
damping, cf. Ref. |20].) This refiects the fact that out of 
equilibrium, the electronic degrees of freedom can effec- 
tively pump energy into the collective modes. 

(iii) If the conductor is time-reversal symmetric such 
that S — S'^ and A — —A^, the effective orbital mag- 
netic field 7° vanishes in thermal equilibrium. In general 
out-of-equilibrum situations, 7° is nonzero even for time- 
reversal symmetric conductors since the current provides 
another quantity which is odd under time reversal. 

(iv) In thermal equilibrium, the variance D,^^' of the 
fluctuating force is related to the dissipation matrix 
7" through the fluctuation-dissipation theorem D,j^' = 
2Tj^^,. This follows readily using that the second term 
in Eq. ( |12| vanishes in equilibrium and that —d^f = 
/(I — f)/T. We note that a recent study of back- 
action forces in quantum measurements in terms of the 
5-matrix obtains a friction coefficient which violates the 
fluctuation-dissipation theorem unless the S-matrix is en- 
ergy independent. 

(v) The current flowing in lead a is given in terms 
of an appropriate average over the Langevin process of 
the current to linear order in X. Starting with the GF 
expression, this is the sum of the conventional Landauer- 
Biittiker current P^, /q(X) = {e/h)^pj de{fa — 
fp) Tr(iS'n^5^nQ), and a pumping correction |18] . 



%:„.T.(n.gn,s. 



+fp ReTr ( *n„— H^— - 2n„An^St 



(15) 



In thermal equilibrium, the first term reproduces the 
well-know S-matrix expression |21| for the pumping cur- 
rent, while the second term vanishes by unitarity of S . 

Limit-cycle dynamics. — We now employ our general 
theory to show that at finite bias, the current-induced 
forces can change the collective dynamics qualitatively. 
Without applied bias, the collective dynamics exhibits 
oscillations about a stable equilibrium. These oscilla- 
tions can become destabilized under out-of-quilibrium 
conditions, due to negative eigenvalues of the dissipa- 
tion matrix 7*, a non-conservative force field F(X), or 
both. Within a stability analysis for these scenarios, one 
would predict an exponential instability of at least one 
collective mode jl5| . 

A more thorough understanding of the resulting dy- 
namics requires one to retain the full anharmonic 
Langevin dynamics, including the complete dependence 
of the current-induced forces on the mode coordinates 
X, which is readily possible within our formalism. In- 
deed, we find that the current-induced instability typi- 
cally drives the collective degrees of freedom into a limit 
cycle. Systems with more than one collective degree of 
freedom are most susceptible to this limit-cycle scenario 
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Figure 1: (Color online) (a) Limit cycle (blue solid line) and its approach (blue dotted line) at large bias vs stable oscillations 
at low bias (red asterisk) in the Langevin dynamics without fluctuating force, (b) Several periods of typical trajectories with 
fluctuating forces for the parameters of the limit cycle in (a), (c) Fourier transform of the correlation function {Xi{t)X2{t + T)) . 
The limit cycle is signaled by a single peak, as opposed to two peaks in the absence of a limit cycle, (d) The same signature 
appears in the current-current correlation function, making the onset of limit-cycle dynamics observable in experiment. These 
zero-temperature results, based on analytical results for the current-induced forces, are obtained for Ai — Xiao and A2 = X^/^x 
with A1/A2 = 3/2. In units of X^/IMujq) [where A — (Ai -I- A2)/2], the elastic modes are degenerate with hujo = 0.014, 
Fl.h = ^^P^ (ao ± (Tz), and the hopping between the orbitals is w; = 0.9. The dimensionless coordinates are Xi — {Mujg/ X)Xi. 



when the elastic contributions to the mode frequencies 
are (near) degenerate. 

Consider a model, inspired by transport measurements 
on hydrogen molecules between electrodes [55], with two 
electronic orbitals which are each coupled to one of the 
leads and which allow tunneling between them. The hy- 
drogen molecule in a single-molecule contact has two al- 
most degenerate low-energy vibrational modes [22', cor- 
responding to the center-of-mass vibration between the 
contacts and the rigid-rotation mode. Taking their re- 
spective couplings Ai^ to the electronic orbitals as Pauli 
matrices in the electronic subspace, say ctq and cr^, we 
find that the mean force becomes non-conservative in the 
presence of a finite bias. Without fluctuating force, this 
gives rise to limit-cycle dynamics at larger biases which 
evolves from the stable equilibrium at small voltages, see 
Fig. 1(a). With Langevin force [Fig. l(b)[, the trajectory 
fluctuates about the limit cycle, with occasional larger 
excursions. 

Most importantly, this transition into limit-cycle dy- 
namics has directly observable consequences even for the 
full Langevin dynamics. In the absence of a limit cy- 
cle, the collective motion is a random superposition of 
two normal modes whose frequencies differ due to the 
current-induced forces. (Here, we assume for simplic- 
ity that the elastic contributions to the mode frequen- 
cies are strictly degenerate.) Apart from higher harmon- 
ics, this leads to two distinct peaks in the Fourier trans- 
forms of correlation functions such as {Xi{t)X2{t + r)). 
In contrast, once the limit cycle develops, its frequency 
completely dominates the correlation function, and its 
Fourier transform exhibits a single peak only [Fig. l(c)[. 
This frequency locking is mirrored in the correlation func- 
tion of the current through the NEMS, as shown in Fig. 
1(d). This should serve as a direct signature of current- 
induced limit-cycle dynamics in experiment. Remark- 



ably, the present system resembles a nano-scale motor as 
the mechanical motion remains approximately periodic 
even for the full Langevin dynamics and this periodic 
motion is driven by a non-conservative force. 

We acknowledge discussions with P. Brouwer, G. 
Zarand, and L. Arrachea as well as support by the DFG 
through SPP 1459, SFB TR/12 and SFB 658. 
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